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Summary 

The  imaging  system  conventionally  used  for  Com¬ 
puterized  Ionospheric  Tomography  (CIT)  [1]  is  an  ex¬ 
ample  of  a  limited  angle  tomography  system.  The  aim 
is  to  image  the  electron  density  distribution  of  the  iono¬ 
sphere  in  the  two-dimensional  plane  lying  between  a 
satellite  orbit  at  an  altitude  of  1000  km  and  a  co-planar 
chain  of  ground  station  receivers  as  shown  in  Figure  1 . 


The  ionosphere  lies  approximately  in  the  80  —  500  km 
altitude  range.  As  the  satellite  passes  over  head  it  sends 
dual  frequency  radio  signals  at  150  and  400  MHz  re¬ 
spectively  which  undergo  different  amounts  of  Doppler 
phase  shift  along  each  ray  path  before  they  are  simulta¬ 
neously  received  by  the  ground  stations.  The  difference 
of  Doppler  shift  encountered  by  each  received  signal 
sample  is  subsequently  processed  to  obtain  informa¬ 
tion  about  the  total  electron  content  (TEC)  of  the  iono¬ 
sphere  along  the  corresponding  ray  path.  TEC  records 
for  parallel  ray  paths  are  grouped  together  as  samples 
of  pseudo  tomographic  projections,  and  used  to  recon¬ 
struct  the  ionospheric  plane  electron  distribution  with 
the  aid  of  different  tomographic  image  reconstruction 
algorithms  [3,4,5].  The  result  of  missing  horizontal 
rays  is  that  reconstructions  lack  considerable  distribu¬ 
tion  information  in  altitudinal  direction  [6].  Hence,  it  is 
difficult  to  locate  absolute  heights  of  important  electron 
formations  such  as  peaks  and  troughs  [2]  as  shown  in 
Figure  2. 
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Figure  1.  CIT  imaging  system  geometry. 
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Figure  2.  Image  reconstruction  without  and  with  volumetric  correction. 
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Many  current  reconstruction  algorithms  use  infor¬ 
mation  from  model  ionospheres  such  as  IRI-90  [2]  to 
supplement  the  measured  TEC  data  for  reconstructions 
with  better  shape  definition  and  background  support, 
while  several  others  use  a  priori  free  reconstructions 
involving  smoothness  constraints  and  extremal  nullity. 
Many  of  these  algorithms  employ  a  priori  or  mathe¬ 
matical  constraints  on  a  global  basis  with  little  scope 
for  regional  adjustments  or  allowances  for  unique  and 
unexpected  features  of  the  measured  ionosphere  which 
are  not  completely  representable  with  known  structural 
or  behavioral  information. 

The  Volumetric  Constraints  approach  to  CIT  im¬ 
age  reconstruction  introduces  weighted  correction  val¬ 
ues,  based  on  local  deviations  of  the  current  image  from 
a  priori  parameters,  during  each  cycle  of  a  simultane¬ 
ous  iterative  reconstruction  (SIRT)  process.  The  funda¬ 
mental  and  modified  steps  of  Volumetric  SIRT  (VSIRT) 
reconstruction  are  shown  in  Figure  3.  Volumetric  or 
shape  correction  of  the  reconstruction  is  performed  at 
two  stages  during  a  given  iteration  of  the  algorithm. 
The  first  instance  occurs  during  the  redistribution  of 
the  difference  between  TEC  values  and  those  calcu¬ 
lated  for  corresponding  rays  from  the  currently  known 
image.  The  second  instance  involves  volumetric  shap¬ 
ing  in  the  image  domain  after  TEC  based  correction  has 
been  made.  The  shape  of  an  image  is  defined  to  be  the 


normalized  convolution  product  between  the  image  and 
some  form  of  smoothing  mask.  Ionospheric  electron 
distribution  patterns  typically  yield  very  smooth  low 
frequency  images,  hence  the  mask  chosen  for  shape 
extraction  is  a  large  circular  smoothing  template  which 
preserves  large  scale  significant  variations  while  elim¬ 
inating  minor  local  differences  that  might  arise  in  in¬ 
dividual  images. 

This  corrective  action  in  both  the  projection  and 
image  domains  in  each  iteration  aims  to  achieve  a 
balance  between  one-dimensional  distribution  patterns 
along  individual  ray  paths  and  two  dimensional  ones 
suggested  by  the  position,  extent  and  electron  density 
of  various  ionospheric  layers.  In  each  instance  a  nor¬ 
malized  weight  vector  is  used  to  simultaneously  rein¬ 
force  the  natural  shape  of  the  reconstruction  and  direct 
it  towards  that  of  the  most  probable  combination  of  a 
priori  features.  Thus  localized  reconstruction  of  fea¬ 
tures  most  favored  by  actual  measurements  is  substan¬ 
tiated  with  relevant  a  priori  support  both  in  the  projec¬ 
tion  and  image  domains.  The  actual  location,  extent 
and  electron  density  of  ionospheric  formations  in  the 
reconstructed  image  is  validated  by  rescaling  relevant 
portions  of  the  image  to  satisfy  simulated  “on  site” 
ionosonde  and  ionospheric  backscatter  sounder  infor¬ 
mation. 

TEC  data  obtained  from  simulated  IRI-90  model 


|  Measured  TEC 
j  samples 


Correction 


Image  j 

< 

r  1 

i  Redistribute  i 

—  ^ 

T 

Image 

estimate  Calculated  TEC 

|  error  along  > 

estimate 

for  current  \  samples 

1  rays  , 

for  next 

iteration  j 

TEC  sample  error 

iteration 

4 


ii 


A  priori  TEC 
samples 


-►< 


1 


Error  based  j 
weight  | 
evaluation  I 


ZEZL 


Profile  based  j 
error  \ 

redistribution  i 

— i - 1 


Correction  Image 

■W+hH 


Corrected] 

Image 


Image 

shape 


TEC 

profiles 


Weighted 

Average 

profiles 


A 


Shape  error 
istribution  in 
pixel  nhbd 


A  priori  Data  Bank 
of  Images 


(  +  ) 

jf +'* 

Pixel 

wise 

comparision 

weighting 

correction 


A  priori 
shapes 


f 


Original  SIRT 
Volumetric  SIRT 


I  Image 
|  estimate 
j  for  next 
I  iteration 


Figure  3.  The  Volumetric  Simultaneous  Iterative  Reconstruction  Algorithm. 
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images  were  used  to  form  reconstructions  based  on  sev¬ 
eral  different  types  of  test  suites.  Figure  4  shows  two 
such  test  suites  and  the  reconstructions  derived  from 
them.  In  each  case,  the  contours  of  the  original  image 
are  shown  in  the  background  as  dashed  lines,  enabling 
comparison  of  reconstructions.  Also,  the  test  image 
contours  show  how  different  the  a  priori  sets  are  from 
each  other  and  the  original,  thus  supporting  the  claim 
that  volumetric  reconstruction,  by  nature,  provides  sub¬ 
stantial  data  directed  conformity.  In  particular  compar- 
ision  of  the  reconstructions  of  Figure  4  with  the  non- 
volumetric  reconstruction  of  Figure  2  shows  the  im¬ 
provement  contour  conformity  achieved  in  the  latter. 
The  vertical  position  and  extent  of  the  peak  between 
20-30°  latitude  is  clearly  defined  in  the  volumetric  re¬ 
construction.  Within  figure  4,  the  a  priori  set  on  the 
left  hand  side  has  more  similarities  with  the  test  im¬ 
age  than  the  right  hand  one.  However,  reconstructions 
relying  on  either  set  show  substantial  conformity  with 
the  test  image,  showing  how  VSIRT  selects  the  most 
likely  volumetric  distribution  pattern  even  from  a  di¬ 
verse  and  dissmilar  suite  of  a  priori  images.  More  test 
cases  will  be  analyzed  and  presented  in  the  complete 
version  of  this  paper. 
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Figure  4.  A  priori  sets  1A  and  1C  with  reconstructions. 
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